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Abstract 

The relativistic Fokker Planck equation has been used to study the evolution of 
the quark distribution in the quark gluon phase expected to be formed in ultra- 
relativistic heavy ion collisions. The effect of thermal masses for quarks and gluons 
is incorporated to take account of the in-medium properties. We hnd that the kinetic 
equilibrium is achieved before the system reaches the critical temperature of quark 
hadron phase transition. We hnd that chemical equilibrium is not achieved during 
this time. We have evaluated the electromagnetic probes of quark gluon plasma 
from the non-equilibrated quark gluon phase and compared them with those in 
completely equilibrated scenario. The hard QCD production rates for the electro¬ 
magnetic ejectiles as well as the heavy quark production rates are also calculated. 

PACS: 25.75.+r,12.38.Mh,24.60.Ky,24.85.+p 

Keywords: Quark-Gluon Plasma, photons, dileptons, diphotons. Brownian mo¬ 
tion. 


I. Introduction 


The theory of strong interactions - Quantum Chromodynamics (QCD) predicts 
that under extreme conditions of large baryon density and/or high temperature, 
the hadronic system would dissolve into their fundamental constituents, the quarks 
and gluons. It is expected that the temperature and density achievable in nucleus- 
nucleus collisions at ultra-relativistic energies would be favourable for the formation 
of such a phase, called the Quark-Gluon-Plasma (QGP). It is of fundamental impor¬ 
tance to understand whether thermodynamic equilibrium is achieved in the quark 
gluon system, so as to justify the concept of “phase”. Ideally, the transport theory. 
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or the kinetic theory, should provide an appropriate framework to consider systems 
out of equilibrium. The application of Boltzmann equation to relativistic quantum 
systems is however laced with major difficulties!]^. The non-abelian nature of QCD 
has made this problem rather formidable and as a result, the ’’transport equation” 
approach to many body QCD has found only limited success to date. Nonetheless, 
the held is in rapid progress and some important results have already been obtained. 

The primary motivation to study the non-equilibrium evolution of the quark- 
gluon system is driven from the fact that the characteristic time scales for the 
partonic processes (g, q and g) are of the order of the lifetime of the putative QGP. 
Even if the system achieves thermodynamic equilibrium at some point of time, the 
study of the pre-equilibrium aspects is important to evaluate in the sense that the 
pollutants from this era may affect the kinematical domains where one looks for the 
signals of QGP. QGP diagnostics rely quite heavily on the phase space densities and 
distributions of quarks and gluons. To what extent equilibrium is achieved should 
obviously affect these signals. 

To this end, the mechanisms governing the approach to thermalisation in the 
quark-gluon system have been a very topical issue of late [0, ^ . Recently a parton 
cascade picture which purports to study the QGD-based space-time evolution of the 
partonic system has received a fair amount of attention!^. The importance of the 
microscopic approach embodied in such pictures notwithstanding, their applicability 
to actual QGP diagnostics is still largely unexplored or beset with numerical diffi¬ 
culties. Also, the approximations inherent in the parton-cascade model have been 
questioned by some authors [^. In this work, we thus propose to use a physically 
transparent, semi-classical model to understand the evolution of the many body 
quark-gluon system towards equilibrium. 

The central theme of our approach is to exploit the well-known result that gg 
cross-section is considerably larger than qg or qq cross-sections, primarily because 
of the colour factor of gluons. It is therefore reasonable to expect that the gluons 
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would thermalise among themselves appreciably earlier than the whole system of 
quarks , antiquarks, and gluons. The proper time Tg at which the gluons equilibrate 
is thus considerably less]^] than the overall equilibration time Tq; the value of Tq was 
proposed to be of the order of 1 fm/c by Bjorken some time ago. 

The gluons carry about 50% of the momentum and sea quarks only a tiny frac¬ 
tion. Thus, in very high energy collisions (RHIC or LHC energies), if we confine 
our attention to the central rapidity region, it is quite natural that from Tg onwards, 
the equilibrated gluons may provide a thermal heat bath for the sea quarks (anti¬ 
quarks). This picture is further justihed by the fact that the sea quark (antiquark) 
density is very low compared to that of the gluons in this region|p]. Thus we are left 
with a system where a relatively small mixture of non-equilibrium degrees of free¬ 
dom (quarks and antiquarks) interact with some equilibrated degree of freedom (the 
‘gluonic’ bath); such processes are known to give rise to Brownian motion^] which is 
governed by the Fokker-Planck (FP) equation[^. QCD being asymptotically free, 
hard collisions involving large momentum transfers are suppressed compared to soft 
interactions and in our picture, thermalisation in the quark-gluon system proceeds 
through many such soft collisions. The FP equation describes, semi-classically, the 
evolution of the many body quark-gluon system in a kinetic theory framework. The 
system under consideration is highly relativistic and presumably at high tempera¬ 
tures. Therefore, account of production and annihilation of qq pairs in the gluonic 
heat bath must be taken. 

The goal of this paper is two-fold. First, we develop the framework of the 
relativistic FP equation to describe the temporal evolution of the quark distribution 
function and then use this information to estimate the effect of the pre-equilibrium 
era on the QGP diagnostics. For this latter purpose, we conhne our attention to 
the electromagnetic probes - photons, photon pairs, and dileptons. It is well known 
that these probes constitute an especially clean class of QGP signals]^; by the 
very nature of their interaction; these probes tend to leave the system of strongly 
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interacting matter with minimal distortion of their energy/momentum. They thus 
carry the information from within the reaction zone rather more effectively, not 
being masked by the details of the evolution process. Obviously, such is not the 


case for hadronic probes which lose the initial information because of their strong 


coupling to the rest of the system. 

The structure of the paper is as follows. In Section II, we derive the FP equa¬ 
tion in a form which is appropriate for our present purpose. In Section III, we 
discuss the production of photons, diphotons, and dileptons from non-equilibrated 
scenario. Section IV is devoted to the discussion of our results and we conclude with 
a summary in Section V. 

We reiterate that our analysis is restricted to the situation where the central 
rapidity region is free of baryon number (complete transparency), a situation which 
is expected to be achieved in RHIC or LHC energies. In these circumstances, the 
collective evolution of the system is governed by the scaling hydrodynamics d la 
Bjorken|0] which we tacitly assume to be the case. 


II. The Fokker-Planck (FP) Equation 


The Boltzmann equation in the relativistically covariant form can be written as 


P^dJ{x,p) = C{f} 


( 1 ) 


where p is the four momentum of the test quark and / is its phase space density. 
C{f} is the collision term. The left hand side of eq. (1) is known as “streaming 
term”. In the spirit of boost-invariance incorporated in the scaling hydrodynamics 
of Bjorken, we assume that the phase space density of the quark is independent of 
X, i.e. the plasma is uniform. Under these conditions the Boltzmann equation reads 


df C{f} 


dt E 



coll 


( 2 ) 
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We can seperate the elastic from the inelastic collision term as follows 



coll 



inel 


coll 


( 3 ) 


First we consider the elastic collisions of the test quark with other quarks, anti¬ 
quarks, and gluons present in the system. For the system under study, the gluons 
are in complete equilibrium, with a density larger than the non-thermal densities of 
the quarks and antiquarks present in the system. The rate of collisions w{p,q) is 
given by 

w{p,(l)= ( 4 ) 

j=<l,Q,9 

where denotes the rate of collisions of a test quark q with other participant parton 
j, i.e. for the reaction jq — jq. In our case the term of eq. (4) dominates over 
the other two terms, because of the paucity of quarks and anti quarks in the system. 
The expression for can be written as 


w^{p,q) = -fj 




( 5 ) 


where is the statistical degeneracy (2x8 for gluons) and fj{q) is the phase space 
density for the species j; Vrei is the relative velocity between the test quark q and 
the other participating parton j and is the relevant cross section for the elastic 
collision. 

We assume that the transition takes place between two states having momenta, 
say, p! and p, where pi — p is very small. In terms of collision rates this means that 
the function w{p,pf) is sharply peaked at p ~ pi. The right hand side of eq. (2) can 
be written as 

el 

= f d^k[w{p + k,k)f{p + k) — w{p,k)f{p)] (6) 

coll 

Expanding w{p -|- k, k)f{p -|- k) by Taylor series we get 

w{j) + k, k)f{p + k)^ w{p, k)f{p) + k ■ ^{wf) + 



5 



Substituting eq. (7) in eq. (6) we get, 


dt 


el 


= / d-^k 


coll 


r ^ 1 


i.wf) 


( 8 ) 


For a one dimensional problem, we can write the above equation as 

el 


dt 




where 


^ [/^2(P)/] 

( 9 ) 


(10) 


(F) is the average force acting on the test particle, and 


h 2 (p) = ^ / d^kw{p, k)k‘^ 
Combining eq. (2) and eq. (9) we get 


{{h?) 

6t 


( 11 ) 


f = |lft(p)/l + |5l/..W/l 


( 12 ) 


This is the celebrated Landau kinetic equation [^, a nonlinear integro - differential 
equation. The nonlinearity is caused due to the appearence of / in /xi(p) and /i 2 (p) 
through w{p, k). It arises from the simple fact that we are studying a collision process 
which involves two particles - it should, therefore, depend on the states of the two 
participating particles in the collision process and hence on the product of the two 
distribution functions. As is evident from the derivation, the equation is valid for a 
weakly coupled system, where the average kinetic energy is large compared to the 
two particle interaction energy. 

We can attain considerable simplicity if we replace the distribution functions 
of the collision partners of the test particle by their time independent equilibrium 
Fermi-Dirac or Bose-Einstein distributions (depending on the statistical nature) in 
eqs. (10) and (11). Then eq. (12) reduces to a linear partial differential equation - 
usually referred to as the Fokker-Planck equation]^. 
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To relate /ii and /i 2 with dynamical parameters relevant for the system nn- 
der stndy, let ns consider the classical eqnation of motion of a particle execnting 
Brownian motion in a heat bath, 

I = ( 1 =*) 

where F(t) = Frit) + F^it). Here Frit) is the rapidly flnctnating part and F^it) is 
the viscons force. Taking the average of the above eqnation and assnming that the 
average of Fr{t) over a long interval of time vanishes, i.e. 


(FAt)} = 0 

{Fdit)) ^ -apV ( 14 ) 


where Op is the dynamical friction parameter containing the dynamics of the problem 
nnder stndy (see next section for details), we can, in the relativistic case (n = 
p!\/p^ + m ^), write 

/ii(p) = -apV = (15) 

VP + 

Assnming that the conpling between the Brownian particle and the bath is weak, 
we have 

P2(p) = = 2ap(vp) (16) 

For the nltra-relativistic case n ~ 1, p ~ T, implying 

/i2(p)- 2apT=2DF (17) 


Using eqs. (12), (15) and (17) we get 


dt 


d 


apPf 


+ D, 




(18) 


dp y i/p^rpr ^J 

This is the Fokker Planck eqnation describing the evolntion of a qnark towards 


eqnilibrinm due to its interaction with the gluonic heat bath (see ref. \12 

ollisi 

ay 


The relativistic FP equation with inelastic collisions can be written as 

inel 


A 

dt dpz 


OipPzf 

IpI + rriF p 




( 19 ) 
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We can linearize the above equation with the relaxation time approximation [^, 
n as follows, 

on j - Jeq ^ 20 ) 


f-L 


eq 


dt J Trelax 


where f^g is the equilibrium distribution and Treiax is the relaxation time estimated 
from the reactions gg qq and g <-*• qq, rriT is the transverse mass (= + mljf). 

meff is the effective mass dehned as 


me/f = \J + mLnnoJ 


( 21 ) 


where rricurrent is the current quark mass (= lOMeV for u and d quarks) and rrithermai 
is the thermal mass: 

mthermal = \J 9 sT'^/Q ( 22 ) 

gs is the strong coupling constant. In a chemically non-equilibrated scenario, the 
thermal mass is replaced by ni^hermai = (1 + '’^ql‘^)9lT‘^ [0]’ where Vg is the ratio 

of equilibrium to non-equilibrium density. We have seen |jT^ that the effect of such 
a change on thermal mass has negligible effects on the hnal results. 

The FP equation reduces to. 


A 

dt dpz 


(^pPzf 


D, 




/-/, 


eq 


(23) 

have 


/p'^ -|- J ^Pz T'relax 

It should be mentioned at this point although several authors 
used the relaxation time approximation to study the approach to equilibrium in a 
quark-gluon system, such a treatment is meaningful only for small deviations from 
the equilibrium conhguration. We have dwelt on the relaxation time approach in 
some detail only to clarify the physical picture. A more consistent way is to evaluate 
the contribution of the inelastic term through a time-dependent normalization of /, 
which can be estimated by solving the momentum-integrated Boltzmann equation, 
taking proper account of the reactions g ^ qq and gg <-*• qq. These details 


have been reported in the literature [18|; for the sake of brevity, we do not repeat 
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them here but refer the reader to this work. It should be noted that the reactions 
9 9 9 9 9--- do not appear explicitly as the gluons have been assumed to be 

thermalized so that their density is determined from the temperature of the bath. 
Also, the thermal mass of the gluons is an essential ingredient; otherwise the reaction 
g ^ qq would be forbidden. 


Ila. Determination of a. 


We assume that the phase space distribution function can be factorised as f{p, r) = 
f{Pz,T)G{pT), where r is the proper time, G{pt) = exp(—p|,//i^)/7r/i^ and f{pz,T) 
is the solution of FP equation. 

The friction parameter Op is a very crucial parameter. It contains the dynamics 
of elastic collisions {qq,qq and qg). It can be related to the energy loss of a quark 
in a dense partonic medium ||^ in the following way: 


ap{Pz,r) = 


Pz 


dx I 


ap(r) = ( 


E 

Pz 


dE' 

dx 


) 


(24) 


The approach to equilibrium for the different quark species is then determined by 
eq.(23). In principle, Op may be determined from kinetic theory formulation of 
QCD through the application of the fluctuation dissipation theorem |T0|, but that is 
indeed far too complex a problem to handle, given the present state of the art. It 
can, however, be assumed that since the friction constant is expected to be largely 
determined by the properties of the ’’bath” and not so much by the nature of 
the test particle, one may take ap(p^, r) ~ apij). In this respect, we recall the 
earlier work of Svetitsky |pO|| where the classical diffusion and drag coefficients of a 
nonrelativistic charm quark propagating in a quark gluon plasma were calculated. 
Although his scenario is somewhat different from ours, the operating equation in 
both cases happens to be the Fokker Planck. In his dynamical calculations, he found 
approximate momentum independence of the drag coefficient (fig.2 of |^), entirely 
in line with our assumption. It is however not realistic to the values of Op from 
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for lighter quarks. We may also remark here that a recent work has appeared in 
the literature where a Fokker Planck type equation, including the non-abelian 
features of QCD in the collision terms of the transport equation, has been discussed. 
The main attraction of this work is in studying the damping of the collective colour 
modes, of relevance to jet quenching studies but outside the scope of the present 
work. There is also a component which governs diffusion in momentum space, but 
the deviation from the abelian case is rather small. The correction is proportional 
to the small non-equilibrium deviations and as such can be generally neglected 


It is however noteworthy that these authors also relate the momentum diffusion 
(or friction ) constant to the partonic dE/dx, as in the present work. Let us also 
mention that we have assumed the temperature T(r) to arise from the thermal bath 
whereas these authors look at non-equilibrium contributions to both fg and fq { fq). 
There have nonetheless been some recent developments [^, ^ in connection with 
jet quenching studies in QGP which may shed light on this issue. The expression 


for energy loss has been calculated by various authors in recent times and for 
light quarks, it is given by 




if — 1 , 1 -I- n' 

n H-h In ■ 


(25) 


dx 3 ^ 2 ' ^ ~ 'iJ, 

where ag is the strong coupling constant, Gi? is the Casimir invariant, k^ax is fhe 
maximal momentum (~ p, the momentum of the light quark) and kD is the Debye 
momentum. The value of ag is calculated from the following parameterisation 


ag = 67r/(33 — 2nf)ln{KT/Tc) For heavy quarks with E « Mq/T the expression 
for dE/dx can be written as 


dE 

dx 




1 + 




X In 


1 1 
ET ' 

rrigMQ^ 


v" , 1 -|- n 


2n2 


In 


(26) 
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For E » Mq/T we have 


dE 

dx 


Svra^T^ 


(l + In |^2’"//2(6 +^/)o.920^^^^ (27) 


where Mq is the mass of the heavy quark, B{v) is a smooth function of velocity 
(see Ref.|2^ for details), rrig is the thermal gluon mass, rrig = 5ffT^/3(l + Nf/6). 
In the region of E ^ Mq/T, we have used the eqs. (26) and (27) appropriately. 


as in Ref. [23], It is important to mention here that we have included the radiative 


energy loss |2^ in dE/dx {dE/dx \rad.= 27ia‘lC2T^ (modulo log terms)); however, 
the effect of such processes is rather small. Recently, Baier et al have calculated 
the energy loss due to radiative process (~ ag^JEqjj/lg \n{E/q]^lg)) including the 
rescattering of the radiated gluons which is ignored in is the mean free path 

of the gluons. Weldon has shown that the energy loss is proportional to for 
radiative process as well as for elastic collisions if one assumes Ig ~ a^T. 


Ilb. Cooling of the Gluonic Heat Bath 

The gluonic heat bath is cooling due to expansion and the rate of cooling is deter¬ 
mined by the relativistic hydrodynamics. The bulk properties of the system, e.g. 
the cooling law etc., are governed by the equilibrated degrees of freedom. In our case 
the cooling law is given by the Bjorken’s scaling law with appropriate modihcation 
due to quark production through thermal gluon fusion gg —> qq and thermal gluon 
decay g —>■ qq. We have obtained the cooling law||18|| by solving the hydrodynamic 
equation which is parameterised as T = a/r^ where a = 0.4077, /3 = 0.355 at LHC 
and a = 0.33, f3 = 0.352 at RHIC energies respectively. The cooling rate in Bjorken 
model (/? = 1/3) is slower compared to the present case where the production of 
quarks cost some energy. We would like to mention here that the effects of trans¬ 
verse expansion on cooling law would be negligibly small for the intial parameters 
under consideration for RHIC and LHC energies (see ref.[^ for details). 


11 































lie. Solution of the Fokker-Planck Equation 


We solve the FP equation with the following initial and boundary conditions for a 
quark species j 


fjiPz,T) 


Aj6{pz 


and 


r ( \ IPzHoo 

fj{Pz,T) -^ 0 


(28) 


(29) 


The parameter Aj is determined from the initial density of the partons. 6{pz) is 


a rather good approximation of the low x structure function|^. We should also 
mention here that the hnal outcome of the model is insensitive to the functional 
form of the initial distribution function - a typical characteristic of the Markovian 
process. 

We can either solve the inhomogeneous FP equation (eq. (23)) with the initial 
and boundary conditions given by eqs. (28) and (29), or equivalently, solve the ho¬ 
mogeneous FP equation with a time-dependent normalization for / which accounts 
for the evolution of quark density as a result of the inelastic reactions g qq 
and g g 


qq 


. Our actual calculations show that there is not much difference 


between the two approaches, insofar as the electromagnetic signals are considered. 
We therefore show here the results for the relaxation time approach, as most other 
authors in this held have worked in this framework. 

It is important to mention here that the hnal outcome depends on the value 
of Aj. There is a lack of consensus about the initial value of the quark density. 
We take the initial values of quark densities from HIJING|^, The phase-space 
density of quark is larger in case of parton cascade model [@] and also in the work 
of Shuryak[P]. In this sense our work corresponds to a conservative situation. The 
data from RHIC and LHC should make a distinction among various models. 
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III. Electromagnetic Probes 

Ilia. Single Photons 


As mentioned in the introduction, photons and dileptons are the most efficient sig¬ 
nals for QGP. However, apart from the photons from QGP, there are other sources of 
photons 1 ]^, e.g., photons from hadronic reactions, initial hard collisions of partons 
and hadronic decays. In this section we shall evaluate the photons from a non- 
equilibrated quark gluon system, thermalised QGP, and from initial hard collisions. 

The dominant reaction channel for single photon emission are the annihilation 
{qq gj) and Gompton {q{q)g —>■ q{q)'y)- The transverse momentum distribution 
of photons produced in a reaction (1 -|- 2 —>• 3 -|- 7 ) is given by 


dN 

d’^Prdy 


Af 


16(27r 
X 


;TtR^ 


^ I TdTdripiTdpiTd6idp2Tdyidy2 

\M\‘^ 


\piT sin( 6 'i - 62 ) +Pt sin 02 l^o 

xf{Piz,r)G{piT)f{P 2 z,r)G{P 2 T) (1 ± f{P3z,r)G{p3T)) (30) 


with 


where 


0 ° 


R 

ip 


= Ip 


cos 


-1 


(JL.) 


= tan 


Y 2Rp2T ) 

{p\t + Pt- ‘2ptPit cos 0i)^/^ 
/ Pit sin 01 \ 


.piT COS 01 — Pt 


( 31 ) 


H = m\ + m\ — m\ + 2 pitPt cos 0i -|- 2 miTM 2 T cosh(2/i — P2) 
- 2 miTPT cosh(|/i -y)- 2 m 2 TPT cosh (?/2 - y) 


(32) 


We have written the total distribution function as f{p, r) = G{pT)f{pz, t) where 
G{pt) = exp(—pf.//i^)/ 7 r/r^ with p = 0.42 GeV[^. 
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For the photons from non-equilibrated and eqnilibrated QGP we take the phase 
space distribntion from the solntion of FP eqnation and Fermi-Dirac distribntion, 
respectively, for the qnarks. Glnons are always described by Bose-Einstein distribn¬ 


tion. 


The hard QCD photon spectra from the Gompton process is evalnated by 
nsing the following expressions, 


da^jy = 0) 

dyd^PT 




dx„ 


?>s‘^{xt/2) Jxmin Xa - {xt/2) 
X 


^2{Xa'i ^)G{Xb] B)— - 2 J - 1 " ^ A ^ B) 


x„x 


a-^b 


where 


Xb = 


XaXT 


I ^min — 


XT 


( 33 ) 


( 34 ) 


2Xa - XT 2-XT 

F 2 {,x) = xJ2^q [q{x) + q{x)] and G{x) = xg{x), where g{x) and q(x) (q{x) ) are the 
strnctnre fnnctions for glnons and qnarks (antiqnarks), respectively. 

Similarly, the resnlt for the annihilation graph is given by 

da^{y = 0) Satts A dxa 


dyd^PT 


9s2 


Xa - {xt/2) 


. xi + xi 


Q2{,Xa] A)Q{Xb] B) °' + {Xa ^ Xb] A ^ B) 

ry*0 /y»'J 


( 35 ) 


where we have dehned Q 2 (x) = xj^ e^qi^x), Q{x) = xJ2 q{.x) 
The strong ’’rnnning conpling constant” is given by, 

127r 


asm = 


(36) 


(33 - 2Nf) HQyA^gcn) 
where Nf is the nnmber of flavonrs and Aqqd is the QCD scale parameter. 
We have chosen for evalnation of single photon spectra from hard QCD 

processes. 
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Illb. Dileptons 

The dilepton spectrum has been considered a promising probe for the QGP di¬ 
agnostics. In the dilepton mass window M = 2-5 GeV, the main source is the 
eletromagnatic annihilation of quarks and antiquarks. The invariant dilepton mass 
spectrum is given by 

IKPd-y = 1^. j rird^OPTdprqTdqr 


where T = \po'P — ppo\/■ The basic cross-section for dilepton production due 
to qq annihilation is given by 


cr, 




/ _ 4my 

9 M^\ M2 ; M2; 


X 


2ml 2m‘i 

1 H- - H- - + 

M2 M2 


4m2m^ 

M4 


- 1/2 


(38) 


We obtain f{pz,T) from the solution of FP equation with appropriate boundary 
conditions as mentioned above. The dilepton production rate in complete thermal 
equilibrium (both kinetic and chemical) is obtained by replacing the quark distri¬ 
bution function by Fermi-Dirac distribution function. 

The Drell-Yan production has been obtained by using the expression given below 


11 . 


da 


47ra2 


Nf 


dM‘^dy 


3sN^M^ 




(39) 


The structure functions qf{x) has been taken from ref. |]33 


IIIc. Diphotons 

For kinematic purposes, it is naturally preferrable and useful to work with a pair 
of particles in the final state which would allow an invariant mass identification 
and thus a filter. To this end, working with dileptons or diphotons is advantageous 
compared to single photons. The disadvantage with diphotons is that it is a 
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process whereas single photon is a aag process in the lowest order; consequently a 


substantial decrease with respect to single photon rates. The possibility of treating 
diphoton as a signal of QGP was considered in |^, . The basic cross-section for 

diphoton production is 


<(A/) = 


X 


4ml 

1 + 


8m 


41 


M2 M^ 


In 


M 2 

2m2 


1 + 


4m2 


211 / 2 ' 


1 w 

M2 




M2 


1 / 2 - 


(40) 


The invariant mass distribution for the diphoton is the same as eq. (37), only the 
dilepton cross-section being replaced by diphoton cross-section. We should however 
note that the temperature dependence of the two cross-sections through the thermal 
mass is quite different. For diphoton production from a thermalised QGP and also 
from the initial hard QGD collisions, the calculation proceeds along the same line as 
for dilepton production. For the evaluation of dileptons from DY process and hard 
diphoton we have chosen to be equal to M^. 


IV Results 

IVa. Approach to Equilibrium 

As mentioned earlier, in the system under study the quark density changes with 
proper time due to two mechanisms. The expansion dynamics (flow) dilutes the 
density and on the other hand, the creation of quarks in the relativistic heat bath 
enhances the quark density. The gluon density decreases due to expansion only. We 
calculate the ratio of the width of the distribution in non-equilibrium and equilib¬ 
rium situations, i.e. j The advantage of calculating the ratio is that 

the expansion effects will get cancelled to some extent, though the cooling in the 
equilibrium and non-equilibrium scenario is different as has been mentioned before. 


16 





















It is obvious that the cooling in the non-equilibrium situation should be faster 
than Bjorken boost invariant cooling law: T ~ The reason behind the faster 

cooling in the present scenario is due to the production of quarks in the gluonic heat 
bath through the reaction described above, soaking away energy and thus, cooling 
rapidly. 

In £g. 1 we plot these ratios as a function of the proper time r. At RHIC 
(LHC) the initial thermalisation time, Tg for the gluons is 0.3 fm/c (0.25 fm/c), the 
temperature Tg{Tg) is 500 MeV (660 MeV) and the initial quark density Uu/d is 0.7 
fm“^ (2.8 fm“^). At RHIC energy (fig. 1(a)) we observe that the ratio D saturates 
to a value ~ 1, at a proper time r ~ 3 fm/c, well before the temperature of the 
system reaches to Tc (~ 160 MeV). In fig. 1(b) we plot D for LHC energies; the 
thermal equilibration is complete within the proper time ~ 2 fm/c. Fig. 2 ((a) and 
(b)) brings out the same information in greater detail. We plot here the width of the 
momentum distribution of quarks as well as gluons. The decrease of the width for 
gluons corresponds to the cooling due to expansion. One can readily see the width 
for quarks first increases and then starts falling just like the gluons, indicating a 
clear onset of thermal (kinetic) equilibrium. 

We evaluate the density of quarks in the non-equilibrium scenario by integrating 
the distribution function f{Pz,T)G{pT) over its momentum. The non-equilibrium 
density Uq has an explicit dependence on r and an implicit dependence on r through 
T(r). But the equilibrium density neq{T) has only an implicit dependence on r 
through T(r). The ratio = nipr)/UeqiT{r)) thus assumes an universal feature, 
since the implicit time dependence gets eliminated. The time dependence of the ratio 
Tq can then be used as a ready marker for chemical equilibrium; the time at which 
the explicit time dependence of vanishes, simultaneously with > 1, corresponds 
to the time for chemical equilibration for the flavour q. We observe from fig. (3) that 
Tq neither saturates nor approaches the value unity before the temperature of the 
system reaches the value (160 MeV). Therefore, we conclude that the chemical 
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equlibrium is not achieved in the quark gluon system, although thermal(kinetic) 
equilibrium is. To show the sensitivity of the evolution of on the initial quark 
density we include, in £g. (3), the result from one of our previous calculations |]^ 
obtained by taking initial quark density from the structure functions p, |^. For the 
sake of completeness we also show the cooling law in £g. (3), where qq production 
has been taken into account. We also see in £g. (4) that the production of ss and 
cc pairs from qq [q = u/d), gg fusion and g decay. As expected the results are not 
drastically altered from the case with qq and gg fusion. 


IVb. Electromagnetic Signals 

In the preequilibrium era, the temperature is very high but the quark density is far 
below its fully equilibrated value. In the standard scenario where one assumes the 
system is fully equilibrated at a proper time tq (~ 1 fm/c), the temperature is small 
for a given multiplicity, but the quark density is large. Within the framework of our 
model, we study the effect of these competing aspects on the photon, dilepton, and 
diphoton spectra. Clearly the above effect will be maximum for the reactions which 
involve quark-antiquark annihilation i.e. on dilepton and diphoton spectra. 

In £g. 5 we plot the single photon spectrum for LHC energies. The hard QCD 
photons dominates the spectra for px > 5 GeV, but these photons are under con¬ 
trol through perturbative QCD (pQCD). Photons from non-equilibrated and equili¬ 
brated QGP are not distinguishable upto a of 5 GeV. Low px photons originate 
from the late stage of the evolution, when the temperature approaches the critical 
value Tc = 160 MeV. For px > 7 GeV, the photons from the non-equilibrated QGP 
system dominate, essentially because of the very high initial temperature though 
the system is far from chemical equilibrium. 

At RHIG,( £g .6), the picture is different. Photons from equilibrated QGP 
dominate the spectra upto px = 8 GeV. For px > 8 GeV pre-equilibrium photons 
dominate because of the high initial temperature. At RHIG energies the system 
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remains far from chemical equilibrium, he., the number of quarks and antiquarks is 
rarer compared to their equilibrated value till the critical temperature Tc- So the 
photons are less in number in the pre-equilibrium scenario compared to the equilib¬ 
rium one. The hard photons dominate for pt above 5 GeV, as in the case of LHC 
energies. Once again they are under control through pQCD. We have conhned our 
attention entirely to the annihilation and Compton channels for photon production 
in the quark gluon sector. Even though the number of quarks (antiquarks) is low, 
these channels still dominate over gg g'j, as seen from fig. (7). 

Figs. (8) and (9) show the dilepton count rates at LHC and RHIC energies 
respectively. It is readily seen that at LHC (hg. (8)), preequilibrium dominates over 
equilibrium conhgurations for all invariant masses. For RHIC (hg. (9)), however, 
preequilibrium dominates for M < 2.5 — 3.0 GeV and equilibrium for M > 3.0 
GeV. This can be understood in the following way. At LHC, the departure from 
equilibrium is not so large while at RHIC, it is substantial (hg. (3)). Thus at LHC, 
the higher initial temperature, together with not too low quark density, results in 
preequilibrium emission dominating over the equilibrium scenario for all values of M. 
For RHIC energies, however, the higher initial temperature is largely compensated 
by the very low quark densities. In both cases, Drell-Yan (hard) is seen to dominate 
over both preequilibrium and equilibrium emissions for M >2 GeV but once again, 
they can be taken care of through pQCD. 

Figs. (10) and (11) show the diphoton count rates at LHC and RHIC, respec¬ 
tively. They behave qualitatively in the same way as dileptons. This is to be ex¬ 
pected, since both dileptons and diphotons come from the qq channels and as such, 
they only differ in the elementary cross-sections (including the effects of thermal 
masses). 
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V. Summary and Conclusions 

We have analysed the approach to thermal and chemical eqnilibrinm in a qnark glnon 
system within the framework of a semi-classical, physically transparent model. A 
fnndamental consequence of this picture is that while thermal (kinetic) equilibrium 
is probable, chemical equilibrium is not, even for LHC energies. Even the kinetic 
equilibrium is achieved through a succession of time scales. The central issue of 
this work is to explore the possibility of testing the equilibrium or pre-equilibrium 
scenario through the so-called signals of quark - gluon plasma. 

We have therefore calculated the various electromagnetic probes for such a suc¬ 
cessive equilibration picture. It is seen that both at RHIC and LHC energies, emis¬ 
sion from the preequilibrium phase does affect the electromagnetic signals of the 
quark-gluon-plasma at the kinematic windows thought to be appropriate for such 
studies. If the goal of heavy ion studies is to look for fully equilibrated quark- 
gluon-plasma, then dileptons and photons at M > 2.5 GeV appear to have a better 
chance of being successful at RHIC. Curiously, even though the system appears to 
approach an equilibrium conhguration much more closely at LHC than at RHIC, the 
dilepton or diphoton signals seem to perform better at RHIC energies, at least for 
the above purpose. It should however be noted that while dilepton and diphoton 
spectra behave very similarly, a single photon spectrum has a markedly different 
behaviour. This obviously is due to the fact that the single photon spectrum has 
contribution also from the Compton {qg gy) channel, in addition to annihilation 
{qq gy) where as dilepton or diphoton are sensitive only to annihilation. Thus 
a correlated measurement of y and yy would shed light on the early evolution of 
the gluon density while simultaneous measurement of and yy would test the 

validity of the mechanism visualised here to study the approach to equilibrium. 

Finally the present calculation dNcc/dy at RHIC energies {Tg = 500 MeV) 
turns out to be ~ 0.2 which compares favourably with that of Shuryak (~ 0.3) p. 
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Inclusion of thermal masses suppresses the heavy quark production in the present 
case compared to that of Shuryak’s calculation where the thermal quark masses 
have been neglected. 


We note that in Ref. the authors observe photons from chemically non- 
equilibrated partonic gas is less than that from equilibrated plasma by a factor of 
10“^(10“^) at RHIC (LHC) energies for 1 < < 3 GeV, whereas Shuryak [|] hnds 

an enhancement in both photon and dilepton productions over the equilibrium sce¬ 


nario at RHIC energies. In Ref. P, the kinetic equilibrium is assumed throughout 
the evolution history unlike the present case. However, in the present case at LHC 
photons from non-equilibrated quark-gluon system and fully equilibrated QGP are 
not distinguishable, whereas at RHIC the equilibrium photon yield is order of mag¬ 
nitude larger than that of non-equilibrated scenario. This is largely because of the 
fact that at LHC the initial partonic gas has more chance to get closer to chemical 
equilibrium compared to RHIC. The difference in photon/dilepton spectra from all 
these calculations could be largely attributed to the initial conditions where much 
work needs to be done [Q. The data from RHIC and LHC will make distinction 
among various initial conditions. 

We gratefully acknowledge helpful discussions with Gerry Brown, Edward 
Shuryak and Larry McLerran at various phases of this work . 
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Figure Captions 

Figure 1. Ratio {D) of the width of the non-equilibrium to the equilibrium momen¬ 
tum distribution for u and d quarks as a function of proper time r (a) at RHIC 
energies, (b) at LHC energies. 

Figure 2. Width of the momentum distributions of quarks {u/d) and gluons as a 
function of proper time r (a) at RHIC energies, (b) at LHC energies. 

Figure 3. The evolution of non-equilibrium density normalised to equlibrium density, 
and temperature as a function of proper time r (a) at RHIC energies, (b) at LHC 
energies. 

Figure 4. Heavy quark production at RHIC and LHC energies as a function of 
temperature at y = 0. 

Figure 5. Transverse momentum distribution of single photons at LHC energies. 
The initial parameters for the equilibrium case are taken as Tj = 290 MeV and 
Tj = 1 fm/c. 

Figure 6. Same as £g. 5 at RHIC energies. The initial parameters for the equilibrium 
case are taken as = 250 MeV and Tj = 1 fm/c. 

Figure 7. Transverse momentum distribution of single photons from the channel 
99 91 compared to qq g'y at RHIC and LHC energies. 

Figure 8. Invariant mass distribution of muon pairs at LHC energies with initial 
parameters same as fig. (5). 

Figure 9. Same as fig. 8 for RHIC energies with initial parameters same as fig. (6). 
Figure 10. Invariant mass distribution of photon pairs at LHC energies with initial 
parameters same as fig. (5). 

Figure 11. Same as fig. 10 for RHIC energies with initial parameters same as fig. (6). 
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